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Abstract 

B — > ixlv decay is studied in the effective theory of heavy quark with infinite 
mass limit. The leading order heavy flavor-spin independent universal wave 
functions which parametrize the relevant matrix elements are evaluated via 
light cone sum rule method in the effective theory. The important quark 
mixing matrix element \V U \,\ is then extracted via B — > ttIu decay mode. 
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I. INTRODUCTION 



Weak decays of charmed and beautiful hadrons are quite favorable in particle physics 
because of their usage in determining fundamental parameters of the standard model and 
testing various theories and models. Among these heavy hadron decays the semileptonic 
decays B — > tiIu and B — > plv have been observed experimentally. These exclusive decays 
provide one of the main channels to determine the important CKM matrix element \V u b\. 

The difficulty in studying B — > ttIv and B — > plv decays mainly concerns the calculation 
of the relevant hadronic matrix elements of weak operators, or, equivalently, the correspond- 
ing form factors which contain nonperturbative contributions as well as perturbative ones 
and are beyond the power of pure QCD perturbation theory. Up to present these form 
factors are usually evaluated from lattice calculations, QCD sum rules and some hadronic 
models. 

Sum rule method has been applied to B —> n(p)lv decay in the full QCD and provided 
reasonable results |I|-£§. Since the meson B contains a single heavy quark, it is expected that 
its exclusive decays into light mesons may also be understood well in the effective theory 
of heavy quark, which explicitly demonstrates the heavy quark spin-flavor symmetry and 
its breaking effects can systematically be evaluated via the power of inverse heavy quark 
mass l/m,Q. The effective theory of heavy quark has been widely applied to heavy hadron 
systems, such as B decays into heavy hadrons via both exclusive and inclusive decay modes. 
There are two different versions of effective theory of heavy quark. One is the heavy quark 
effective theory (HQET), which generally decouples the "quark fields" and "antiquark fields" 
and treats one of them independently. This treatment is only valid when taking the heavy 
quark mass to be infinite. In the real world, mass of quark must be finite, thus one should 
keep in the effective Lagrangian both the effective quark and effective antiquark fields. 
Based on this consideration, a heavy quark effective field theory (HQEFT) j§-[5| has been 
established and investigated with including the effects of the mixing terms between quark 
and antiquark fields. Its applications to the pair annihilation and creation have also been 
studied in the literature [10-12|]. Though the HQEFT explicitly deviate from HQET from 



the next-to-leading order, these two formulations of effective theory trivially coincide with 
each other at the infinite heavy quark mass limit. In our knowledge the exclusive heavy to 
light (pseudoscalar) decay channels have been discussed in [|T^], where the matrix elements 
in the effective theory have been formulated, but the two leading order wave functions have 
not been calculated. 

In this paper we focus on the calculation of the leading order wave functions of B — > nlv 
decay by using the light cone sum rule in the effective theory of heavy quark. As an important 
application, \V u b\ is extracted. In section 2, the heavy to light matrix element is represented 
by two heavy quark independent wave functions A and B. In section 3, we derive the light 
cone sum rules for the calculation of A and B. In section 4, we present the numerical results 
and extract \V u b\. Our short summary is drawed in the last section. 

II. B nLu DECAY MATRIX ELEMENT 

The matrix elements responsible for B — > nlv decay is < Tr(p)\u , y fl b\B >, where b is the 
beautiful quark field in full QCD. It is generally parametrized by two form factors as follows, 
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< 7r(p)\uYb\B(p + q) >= 2/ + (gV + (Utf) + f-tfW- (2-1) 

In the effective theory of heavy quark, matrix elements can be analyzed order by order 
in powers of the inverse of the heavy quark mass 1/ttjq and also be conveniently expressed 
by some heavy spin-flavor indenpendent universal wave functions [p],p],p],[T3|] . 

Here we adopt the following normalization of the matrix elements in full QCD and in 
the effective theory 

— L < n{p)\uTb\B >= -^{< ir(p)\uTQ+\B v > +0{l/m b )}, (2.2) 
where A B = m B — m&, and 

A = lim A B 

m,Q—>oo 

is the heavy flavor independent binding energy reflecting the effects of the light degrees of 
freedom in the heavy hadron. is the effective heavy quark field in effective theory. 
Associate the heavy meson state with the spin wave function 

_ + i I ~7s f° r pseudoscalar meson ^ 
2 1 ( for vector meson with polarization vector e 1 ' 

we can analyze the matrix element in effective theory by carrying out the trace formula : 

< ir{p)\uTQ+\B v >= -Tr^v.p^M,] (2.4) 

with 

ir(v,p) = 7 5 [A(> -p,fi) +pB(v -p,fj,)], 

p» = — . (2.5) 
v ■ p 

A and B are the leading order wave functions characterizing the heavy-to-light-pseudoscalar 
transition matrix elements in the effective theory. They are heavy quark mass independent, 
but are functions of the variable v ■ p and the energy scale [l as well. Nevertheless, since 
the discussion in the present paper is rrestricted within the tree level, we neglect the \x 
dependence from now on. 

Combining eqs. (|2.1|) - (|2.5|) , one gets 



f±(Q 2 ) 



'" bX {A(vp)±B(v P )^} + ---, (2.6) 



\ rn b A 



B V-p 



where the dots denote higher order 1/toq contributions which will not be taken into account 
in the present paper. Note that we have used different variables for / + , /_ and A, B. The 
relation between the variables v ■ p and q 2 is 

m\ + m\ - q 2 

y = vp= . 2.7 

zrriB 
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III. LIGHT CONE SUM RULE FOR B -> -kLu 

The QCD sum rule based on short distance expansion has been proved to be quite 
fruitful in solving a variety of hadron problems. Nevertheless, it is also well known that this 
method meets difficulties in the case of heavy to light transition because the coefficients of 
the subleading quark and quark-gluon condensate with the heavy quark mass terms grow 
faster than the perturbative contribution, which implies the breakdown of the short distance 
operator product expansion (OPE) in the heavy mass limit. Alternatively, it has been found 
that heavy to light decays can be well studied by light cone sum rule approach, in which 
the corresponding correlators are expanded near the light cone in terms of meson wave 
functions. In this way the nonperturbative contributions are embeded in the meson wave 
functions instead of the vacuum condensates in the short distance OPE sum rule. Though 
there are some differences in the techniques of calculation, the two sum rule methods are 
based on the same idea of quark-hadron duality and dispersion relation, and furthermore, 
they follow the same procedure in deriving form factors. 

For B — > irlv decay, one may consider the vacuum-pion correlation function 

F»(p,q) = i J d*xe iqx < Tr{p)\T{u(x)Yb(x),b(0)i^ 5 d(0)}\0 > . (3.1) 

Here p and q are momenta carried by the pion and leptons. The B meson has momentum 
Pb — P + q- Inserting a complete set of states with B meson quantum numbers, we obtain 
the phenomenological representation 

F^(p q) ~ < ^\ U ^ b \ B >K gjgT^O > ] y < 7r(p)\u^b\H >< H\U^d\U > 

m B -(p + q) 2 x m 2 H -(p + q) 2 



With the normalization relation in (|2.2| ), the matrix elements in (|3.2| ) can be expanded into 
the ones in effective theory of heavy quark in powers of l/m&. When all higher l/m* order 
contributions are neglected, (|3.2|) reduces straightforwardly into 



2iF^ — — — ^— + / ds— — — - + Subtractions, (3.3) 
2A B - 2v ■ k J s s - 2v ■ k ' v ; 

where k^ is the heavy hadron's residual momentum, k^ = P£ — m^-iA The first term in 
is a pole contribution obtained by using (|2.4[) together with the parametrization 



< 0\qTQ+\B v >= ^Tr[TM v ] (3.4) 

with F being the leading order decay constant of B meson in effective theory M. The second 
term in (|3.3|) is the higher resonance contributions given in the form of an integral over the 
physical spectral density p(v-p, s). Note that the Lorentz indices of the second term in ( |3.3| ) 
are not written explicitly but embeded in p{y ■ p,s). 

On the other hand, the correlator can be calculated and expressed as the form of an 
integration over the theoretic spectral density p(v ■ p, s) t heory, which equals to p(v ■ p, s) 
under the assumption of quark-hadron duality. Namely, the correlator ( |3.1| ) can be written 
as 
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ds ^ V — S ] + Subtractions. (3.5) 
s — 2v ■ k 

Equating fl3~3[) and ( |3~5D yields 

2iF^ — — — ^— = / ds ^ , + Subtractions. (3.6) 
2A B -2v-k Jo s-2v-k v ; 

So the next step involves the calculation of (|3.1|) in the framework of effective theory 
of heavy quark. Substituting the heavy hadron states and heavy quark fields into the 
effective ones in the effective theory, and then performing the corresponding momentum 
shift Pg — rribV^ = k^, we obtain when neglecting higher 1/mq order corrections 

F"(p,q) = i I d 4 xe i{q - mbv) - x < 7c{p)\Tu{x)YQt(x)>Qt(Q)H 5 d{0)\0 > . (3.7) 



In the light cone sum rule approach, one should contract the heavy quark fields and expand 
the correlator into a series in powers of the twist of light cone pion wave functions. These 
light cone wave functions provide an alternative treatment besides the vacuum condensates. 
They have been discussed in detail in many references |],|2|,|14| . Up to twist 4, the pion wave 
functions relevant to B — > ixlv decay are defined as follows, 

< 7r(p)|iZ(a;)y*7 B d(0)|0 > = -vpf* C due^^u) + x 2 9l (u)] 

Jo 

r 2 fj, r i 



+ - —) / due iup - x g 2 (u) 

x ■ v Jo 



x ■ p Jo 

< 7r(p)|M(x)«7 5 rf(0)|0 > = J*™" m £ due iup - x <f> p (u), 
<n(p)\u(x)a flI/1 5 d(0)\0> = i{ Ptl x v - Pv x^) J * * duc^^u). (3.8) 



Q(m u + m d ) Jo 

<f) n is the leading twist 2 wave function. P and <$> a are twist 3 wave functions, while g\ and 
g 2 are wave functions of twist 4. 

Using the propagator / °° dt5(x — vt) for the contraction of the effective heavy quark 
fields, we get 

F»{y,u) = -4^ dt due-^ u ^{v"[tg 2 (u) - i^<f> p {u) - -^y<t>*(u)] 
2 Jo Jo o 

+ p>*y[-i(f> n - it 2 gi(u) g 2 (u) + -^ a {u)]} (3.9) 

y 6 

2 

with uo = 2v ■ k and u w = -, — ^ — r. 

In order to proceed, we perform a wick rotation of the t axis and then apply the Borel 
transformation B^ to (13.91). The result is 



BPF»(y,u) = -iU f due^ u - l \v»[~g 2 {u) - ^ p {u) + 1-^y^u)] 

Jo 1 61 



+p"y[-&r(u) + ^gi{u) + ^g 2 {u) - -L^^u)}}. (3.10) 



In deriving this equation we have used the feature of Borel transformation: 

4*V» = *(A-±). (3.11) 

It is found that the spectral function used in sum rule can be obtained by performing a 
continuous double Borel transformation on the amplitude itself [|15|,|16|]. In order to get the 
spectral function p(y, s), we now carry out the continuous Borel transformations as follows 



p{y, s) 



M-i/T) £H 



(3.12) 



The result is 

p{y,s) = 



if n Id 

— {^[-^-92^) - /v^>p( 
2y y ou 



u 



1 d 2 
>m) + — — 9i(u) 



y 2 dv? 



?^ 2(M) + %^ (M)]W -*- 



(3.13) 



In the derivation of fl3.13| ), ^ has been first expressed as a derivative of the exponent in 
( P.10| ) over u, and then the method of integration by parts over u has been used. 
( |3lf ) and ( |3.13|) immediately yield: 



m = - 



2A B ~s X Q 

dse^^[——g 2 {u) 
y z ou 



—&Ju) - ^^-0 CT (w)] u= i_^. 
6y du 



2y 



4F Jo 'you y 

B « = -hi" + - 1A*(„) + £A*.(«)U-* , (3.14) 



IV. NUMERICAL RESULTS 



For the light cone wave functions appearing in the sum rules (|3.14 ), we take @,[T^JT7 



V(w) = 6u(l - m){1 + |a 2 [5(2u - l) 2 - 1] + ^-a 4 [2l{2u - l) 4 - U{2u - l) 2 + 1]}, 

2 o 

= 1 + \b 2 \3(2u - l) 2 - 1] + is 4 [35(2w - l) 4 - 30(2w - l) 2 + 3], 
2 8 

3 1 F) 

^(u) = 6u{l - u){l + -C 2 [5{2u - l) 2 - 1] + — C 4 [21(2u - l) 4 - U(2u - l) 2 + 1]}, 

2 o 



gi (u) = -<$V(1 -u) 2 + -e5 2 [u(l -u){2 + 13u(l - u) + 10u d logu(2 - 3u 

+ 10(1 - uf log((l - u)(2 - 3(1 - u) + 2(1 - „)*))], 
g 2 (u)=^5 2 u(l-u)(2u-l). 
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(4.1) 



The asymptotic form of these functions and the scale dependence are given by perturbative 
QCD [0,0. 

For the convenience of comparison, we use the same values for the parameters as in ]2,14 
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a 2 (fi b ) = 0.35, a 4 (// b ) = 0.18, B 2 (fi b ) = 0.29, B^ h ) = 0.58, 
C 2 (ji b ) = 0.059, C 4 (// 6 ) = 0.034, 5 2 (/i b ) = 0.17GeV 2 , e(// 6 ) = 0.36. 



(4.2) 



/ife is the appropriate scale set by the typical virtuality of the beautiful quark, 



m 2 B -m\^ 2.4GeV. (4.3) 



Besides all these parameters the numerical analysis of the sum rules ( |3.14|) needs also 



the hadron quantities /v, f n , A B and F. These quantities have been studied via sum rules 
and other approaches by several groups. With the values /v = 2.02GeV, /„. = 0.132GeV 
A B = 0.53GeV and F = 0.30GeV 3/2 [|, we get from eqs. (|3~T4D the results for A and 



B given in the figures Fig. 1-4. In these figures A and B are shown as functions of T and 
y = v ■ p. We are mainly interested in the range of T = 2.0 ± l.OGeV, where both the twist 
4 corrections and the contributions from excited and continuum states do not exceed 30%. 
It is seen that the curves in Fig.l and Fig. 2 are quite stable in this range for the threshold 
energy s = 2.3 ± 0.6GeV. In Fig. 3 and Fig.4, A and B become rather stable with respect 
to the variation of y = v ■ p when y > 1.5GeV. However, they become unstable at small y, 
which corresponds to large momentum transfer q 2 . This is in expectation because the light 
cone expansion and the sum rule method would break down as q 2 approaches near 0. 

We also derive f+(q 2 ) and f-(q 2 ) from A(v-p) and B(v-p) by using the relations in ( |2.6|) 
and the beautiful quark mass rrif, = tub — A# = 4.75GeV. The results are shown in Fig. 5-6. 
It is readly seen that when the momentum transfer q 2 grows large (e.g. over 16GeV 2 for the 
curve of so = 2.3GeV in Fig. 5-6), the values of /+ and /_ derived from sum rules become 
rather unstable and should not be trusted. 

In order to predict the decay width and \ V u b\, one should have knowledge on the behavior 
of form factors in the whole kinematically accessible region. Now for large momentum 
transfer we have the single pole approximation J3] 

Uiq } = 2m fl .(l- g V<.r ( } 

The couplings fs* and Qb'B-k have been studied in previous papers. Here we would use 
m B * = 5.325GeV, f B * = 0.16 ± 0.03GeV and g B *B-n = 29 ± 3 §. 
Next we write f+(q 2 ) as 

f+(l 2 ) = 1 2 / + 2 ( °l, 4/ 4 ( 4 -5) 

1 — aq z /m B + bq^/m^ 



and fit the parameters a and b by using the sum rules and eq. (|4.4j ). For the threshold 
so = 2.3GeV, we choose proper a and b to make ( [Q| ) approach the sum rule results at 
q 2 < 15GeV 2 but compatible with eq. ( fl4l) at q 2 > 15GeV 2 . Our favorable parameters are 

a = 1.31, 6 = 0.35, /+(0) = 0.35. (4.6) 

The values of f+(q 2 ) at T = 2.0GeV calculated from (fOj), ( |4.5| ) and the light cone 
sum rules are shown in Fig. 7. It is found that the single pole model extrapolation matches 
quite well with the direct estimation from our light cone sum rules ( |3.14 ) at intermediate 
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momentum transfer around q 2 = 15GeV 2 . This implies that our discription of f+(q 2 ) by 
( |4.4| ) together with the sum rules ( |3.14| ) (but in different applicable regions) is self-consistent. 

For the lepton I = e or fi, the lepton mass mi may be safely neglected, and the decay 
width of B — > ttIu has the distribution on momentum transfer q 2 as follows 

J-p ri2 It/ 12 

^ = ^r- (El -™l? /2 \f (4-7) 

Here E n = y = (m 2 B + m 2 — g 2 )/(2me) is the pion energy in the B meson rest frame. 

With the pion mass m n = 0.14GeV and the parametrizations of ( |4.5| ), we obtain the 
integrated width 

T(B -> ixlv) = (10.2 ± 1.5)|Kb| 2 ps -1 - (4-8) 

The error in eq. ( |4.8|) results from the variation of the threshold energy in so = 1.7 — 2.9GeV. 

From the branching fraction measured by CLEO collaboration ||20|| , Br(i?° — > tt~1 + i>i) = 
(1.8±0.4±0.3±0.2)xl0~ 4 and the world average of the B° lifetime |H|,t b o = 1.56±0.06 ps, 
one has || 

T^ -> 7r"/+z/ z ) = (1.15 ± 0.35) x lO^ps -1 . (4.9) 

Comparison of (|4.8|) and ( |4.9|) yields 

|K 6 ] = (3.4 ± 0.5 ± 0.3) x 10" 3 , (4.10) 

where the first (second) error corresponds to the experimental (theoretical) uncertainty. Here 
the theoretical uncertainty is mainly considered from the threshold effects. In general, higher 
order contributions need to be included for all the relevant parameters. It was noticed that 
the two-loop QCD perturbative correction may be significant for an accurate determination 
of B meson decay constants ||. In particular, it may enlarge the constant F by about 
25%, and increase A at the same time. These effects evidently worsen the accuracy of our 
extraction of \V U b\. By taking into account this uncertainty, we arrive at the following result 

\V ub \ = (3.4 ± 0.5 ± 0.5) x 10~ 3 , (4.11) 

This estimate is in good agreement with that derived from full QCD calculation 0: 

\V ub \ = (3.9 ± 0.6 ± 0.6) x 10~ 3 (via B -> nlu), 

\V ub \ = (3.4 ± 0.6 ± 0.5) x 10~ 3 (via B -> plv). (4.12) 
Furthermore, the value of \y u h\ obtained in eq. ( |4.11| ) is also close to the one given by CLEO 



\V ub \ = (3.25 ± 0.14+°;^ ± 0.55) x 10" 3 , (4.13) 

which is a combined result from the analyses based on different models and treatments on 
B — > Tr(p)lu transitions. 
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V. SUMMARY 



In this paper we have studied B — > nlu decay by using the light cone sum rule approach 
within the framework of effective theory for heavy quark. Two leading order wave functions 
in the effective theory with infinite mass limit have been calculated. The important CKM 
matrix element has been extracted and its value has been found to be 

| Kb I = (3.4 ±0.5 ±0.5) x 1(T 3 . (5.1) 

It has been seen that the value of \V u b\ extracted from the leading order heavy quark ex- 
pansion coincides well with that extracted from the full QCD calculation, which shows the 
reliability of the heavy quark expansion and the power of light cone sum rule approach in 
studying heavy to light exclusive decays. Working out 1/mg contributions should be inter- 
esting, and it is expected to cast more light on the treatment of heavy to light decays by 
applying for the effective theory of heavy quark. 
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Fig. 1-2. Variation of A and B with the Borel parameter T for different values of the contin- 
uum threshold sq. The dashed, solid and dotted curves correspond to sq =1.7, 2.3 and 2.9 GeV 
respectively, y = v ■ p = 2.64 GeV is fixed, which corresponds to q 2 = OGeV 2 . 
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Fig. 3-4. A and B as functions of y = v ■ p for different values of the continuum threshold sq. The 
dashed, solid and dotted curves correspond to so =1-7, 2.3 and 2.9 GeV respectively. T=2.0 GeV 
is fixed. 
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Fig. 5-6. Variation of /+ and /_ with respect to the momentum transfer q 2 for different values of 
the continuum threshold so- The dashed, solid and dotted curves correspond to so =1.7, 2.3 and 
2.9 GeV respectively. T=2.0 GeV is fixed. 
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Fig. 7. Variation of / + from different estimations. The dashed curve is calculated from the sum 
rules (3.14) for so = 2.3GeV and T = 2.0GeV. The dotted curve comes from the single pole model 
(|4.4|). And the solid curve is the result of (|4.5|) , which we used to evaluate the integrated decay 
width and \ V u b\. The dashed line and the solid line almost overlap each other at q 2 < 15GeV 2 . 
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